home *** CD-ROM | disk | FTP | other *** search
/ Games of Daze / Infomagic - Games of Daze (Summer 1995) (Disc 1 of 2).iso / x2ftp / msdos / source / snip9503 / rand2.c < prev    next >
Encoding:
C/C++ Source or Header  |  1995-03-14  |  1.3 KB  |  52 lines
  1. /********************************************************************
  2.           The McGill Super-Duper Random Number Generator
  3.              G. Marsaglia, K. Ananthanarayana, N. Paul
  4.   
  5.     Incorporating the Ziggurat method of sampling from decreasing
  6.               or symmetric unimodal density functions.
  7.                       G. Marsaglia, W.W. Tsang
  8.   
  9.                  Rewritten into C by E. Schneider
  10.  *********************************************************************/
  11.   
  12. static unsigned long mcgn, srgn;
  13.   
  14. #define MULT 69069L
  15.   
  16. void rstart (long i1, long i2)
  17. {
  18.       mcgn = (unsigned long)((i1 == 0L) ? 0L : i1 | 1L);
  19.       srgn = (unsigned long)((i2 == 0L) ? 0L : (i2 & 0x7FFL) | 1L);
  20. }
  21.   
  22. long uni(void)
  23. {
  24.       unsigned long r0, r1;
  25.  
  26.       r0 = (srgn >> 15);
  27.       r1 = srgn ^ r0;
  28.       r0 = (r1 << 17);
  29.       srgn = r0 ^ r1;
  30.       mcgn = MULT * mcgn;
  31.       r1 = mcgn ^ srgn;
  32.       return (r1 >> 1);
  33. }
  34.   
  35. long vni(void)
  36. {
  37.       unsigned long r0, r1;
  38.  
  39.       r0 = (srgn >> 15);
  40.       r1 = srgn ^ r0;
  41.       r0 = (r1 << 17);
  42.       srgn = r0 ^ r1;
  43.       mcgn = MULT * mcgn;
  44.       r1 = mcgn ^ srgn;
  45.       return r1;
  46. }
  47.  
  48. /* 
  49. "Anyone who consider arithmetic means of producing random number is,
  50.  of course, in a state of sin" - John Von Neumann
  51. */
  52.